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ABSTRACT 

We test tree-level perturbation theory for Gaussian initial conditions with power 
spectra P(k) oc k n by comparing the probability distribution function (PDF) for the 
density predicted by the Local Lagrangian Approximation (LLA) with the results 
of numerical gravitational clustering simulations. Our results indicate that our 
approximation correctly reproduces the evolved density PDF for — 3 < n < — 1 power 
spectra up to the weakly nonlinear regime, while it shows marginal agreement for 
power indices n = and +1 in the linear regime and poor agreement beyond this 
point. This suggests that tree-level perturbation theory (as realized in the Local 
Lagrangian Approximation) can accurately predict the density distribution function 
for —3 < n < — 1 but fails for n > 0. 

Subject headings: galaxies: clustering, large-scale structure of universe 



1. Introduction 

One of the focal points in the study of large scale structure has been the evolution of P(p), 
the one-point probability distribution (PDF) of the density field (Kofman et al. 1994, Juszkiewicz 
et al. 1994, Bernardeau and Kofman 1995, Protogeros and Scherrer 1996). Whereas in the linear 
regime, assuming Gaussian initial conditions, P{p) scales up self-similar ly by D(t), the growing 
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mode solution, this is not the case once entering the weakly non-linear regime, defined by density 

contrast a = (5 2 = f^^J ) ~ 1- This behavior may be attributed to coupling of different 

Fourier modes, which in the linear regime evolved independently, due to effects of non-local 
interactions in the density field evolution. Furthermore, multistreaming is expected to become 
important at this stage, contributing to the non-linear evolution of the PDF. 

Major progress has been made in the past few years on understanding the quasi-linear 
evolution of the PDF, i.e., the evolution for a <; 1. In particular, a formalism has been developed 
by Bernardeau (1992) to calculate, to lowest order, the hierarchical amplitudes of the evolved 
density field (see the next section for a discussion). More recently, Protogeros and Scherrer (1996, 
PS hereafter) have derived an approximate method (the Local Lagrangian Approximation or 
LLA), which provides a simple analytic expression for the evolved PDF which reproduces, nearly 
exactly, the tree- level hierarchical amplitudes of the "true" evolution. In PS, the predictions of 
the Local Lagrangian Approximation were compared with the results of the "exact" gravitational 
evolution as calculated numerical in a gravitational clustering code. For Gaussian initial conditions 
with tophat smoothing and a scale free n = — 1 initial power spectrum, the LLA predictions 
were found to be in excellent agreement with the "true" evolution. In this paper, we extend this 
comparison to a range of power spectra: — 3 < n < +1. 

Our motivation for undertaking this study is two-fold: i) On a practical level, the Local 
Lagrangian Approximation appears to provide an amazingly simple description of the evolution 
of the PDF for Gaussian initial conditions. We wish to determine if this method is accurate for 
all Gaussian initial conditions, or only for certain power spectra, ii) From a theoretical point 
of view, the accuracy of the Local Lagrangian Approximation can also be considered a test of 
tree-level perturbation theory, since the LLA reproduces, nearly exactly, all of the hierarchical 
amplitudes of the "true" final density field. Previous studies (references) have examined the 
validity of tree-level perturbation theory and the importance of higher-order contributions for a 
few low-order cumulants of the final density field, such as the variance and the skewness. However, 
the LLA provides a method to test all of the hierarchical amplitudes at once. 

We present a short review of the LLA scheme in Section 2 and a description of the N-body 
simulations and our comparative results in Section 3. Our conclusions are presented in Section 4. 
We find that the Local Lagrangian Approximation (and, therefore, tree-level perturbation theory) 
accurately predicts the evolution of the density PDF for initial power spectra with n < — 1, but 
fails for n > 0. 



2. The Local Lagrangian Approximation 
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The Local Lagrangian Approximation is based on the idea that the density at a Lagrangian 
point q at a time t can be approximated as a function only of t and the initial value of p(q): 

pi^t) = f(p(q,t ),t). (1) 

If we define r] = p/p, where p is the mean density, then we choose a mapping of the form: 

The mapping in equation (||) is related to the Zeldovich approximation in the sense that a = 1 
corresponds to the limit of planar collapse (in which the Zeldovich approximation gives an exact 
description of the evolution of the PDF) and a = 3 corresponds to the spherical collapse in the 
Zeldovich approximation. 

In PS, it was shown that the generating function for the hierarchical amplitudes is related 
in a trivial way to the mapping given in equation (|]). In particular, if we take a = 3/2, the 
hierarchical amplitudes for the final density field produced by the mapping in equation @ will be 
almost exactly equal to the true amplitudes produced by exact evolution (see PS and earlier work 
by Bernardeau 1992 and Bernardeau & Kofman 1995). The resulting PDF, P(rj), which is derived 
by applying equation (||) to an initially Gaussian density distribution, automatically satisfies the 
normalization condition / P(rj)ridr] = 1, but it fails to satisfy / P{rj)dr] = 1; this is related, at 
some level, to the problem of multistreaming (see PS for the details). To correct this problem, we 
multiply equation @) by a time-dependent normalization factor N(t): 



N(t) 



leading to the LLA mapping: 



7](q,t) 



f?(q> *). 

<(l-2^(q,i)/3)> 3/2 



|l-25Kq,i)/3| 3/2 ' 
where 5i, the linear-evolved initial density fluctuation, is simply 5i 



(3) 



(4) 



D(t)5 . 



This expression gives the final density for an unsmoothed density field, and so cannot be 
compared directly with either observations or gravitational clustering simulations. If we smooth 
the final density field with a spherical tophat window function, then we obtain a new "smoothed" 
local Lagrangian mapping, f s , given by (Bernardeau 1994; PS): 



Vs = fs(Si) = f 



Sif s (6i 



(n + 3) 
6 



(5) 



where we have assumed a power-law power spectrum P(k) oc k n , and the smoothed mapping 
f s must then be multiplied by the normalization factor given in equation (||). Note that the 
mapping given in equation (||) is not the smoothed density field; rather, it is a density field which 
is guaranteed to give the same hierarchical amplitudes as the "true" evolved density field. 



- 4 - 



Combining equations (Q) and (||), we obtain, for Gaussian initial conditions with spherical 
tophat smoothing, the final PDF: 



(6) 



where the function g is given by 



g{x) 



a 
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Here a = 3/2, o\ is the linear-evolved rms fluctuation: o\ = D(t)ao, and we have defined 

n + 3 



P = 



6 



The normalization factor N is given by: 



N = I g{x)dx. 
h 



(8) 



(9) 



Note that the second term in equation (0) corresponds to the case where the argument of the 
absolute value in equation (||) is negative; this term is negligible compared with the first term 
as long as 01 <J 1. Dropping this second term, we can express the PDF in our Local Lagrangian 
Approximation in the particularly simple form: 



where Si(rj) is given by 



P(rj)dr] 



Si 



1 1 



sf(v) 



N 



N 



and the normalization factor iV is: 

N : 



1 



:x3 N -- 



2<T? 



2-rrai Jo r\ 



(-)e 2 °i d5i{rj). 



(10) 



(11) 



(12) 



Note that equation ( |T2| ) can be integrated by making the change of variables x = tj/N, which 
eliminates N from the integral. 



3. N-body Simulations and Results 
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We calculate the time evolved density PDF in collisionless (Melott et al. 1996) gravitational 
clustering simulations of a pressureless dust, O = 1 universe described more fully elsewhere 
(Melott and Shandarin 1993). We use 128 3 particles on a 128 3 grid and assign initial spectra to 
have the form P(k) = (\5k\ 2 ^ oc k n , with n = —3, —2, —1,0, +1 respectively. The initial spectrum 
high frequency cutoff is given by the Nyquist frequency of the simulation cube ftNy = 64/cj, where 
kf = 2-ir/L is the fundamental mode of the cube and the simulations are terminated at a scale of 
nonlinearity k n \ = 16k f. The initial density contrast variance do is calculated in tophat smoothed 
spheres of radii A = 1,2,4,8, 16 cells and the linear variance corresponding to the final field is 
obtained from 07 = (a n i/ a o) (J Oi where a„;/ao is the ratio of the expansion factors at the end and at 
the beginning of the simulation respectively. We use a cloud-in-cell binning for the final densities 
and consider only 77 in the range [0...4]. The simulation PDF P s (r)) is defined as the fraction of the 
number of volumes of a specific A77 range, Nb, over the total number of grid volumes Nt in the 
simulation, 

ft <"> = i£k- < 13 > 

We expect the rms fluctuation associated with each measured P s (v) value to scale as the 
square root of Nb, provided we take into account only the independent number of volumes arising 
after smoothing, Nj = 3L 3 /47rA 3 , where the L is the side of the simulation box and A is the radius 
of the tophat smoothing sphere. One can then express the rms fluctuation in P s (ri) as: 

ap ^ = iV/Ar/' (14) 



where ctn b = y/NiPjrjjArj . To avoid discreteness effects we only use results corresponding to 
smoothing lengths A > 2 in our comparisons, keeping in mind though that the longer smoothing 
lengths lead to a smaller number of independent volumes within the simulation cube and therefore 
to larger error bars in the evaluation of P s (v)- We examine the evolved PDF results at o\ values 
selected in the range 0.1 < 01 < 2.0 so as to extend our investigation from the linear (07 < 0.5), 
through the weakly nonlinear (07 < 1.0), to the nonlinear (07 > 1.0) regimes. 

Our results are presented in Figure 1, for negative power indices, and Figure 2 for n = 0, +1. 
The solid lines correspond to the Local Lagrangian prediction, while the results of the gravitational 
clustering simulations are given as points with 1 — a error bars. In each graph we quote the 
extrapolated linear variance 07 and the corresponding smoothing scale A. For the n = —3, —2 
models we use a smoothing length 4 < A < 16 whereas for the n = —1, 0, +1 spectral indices we use 
8 < A < 2 in order to examine the PDF behavior in approximately the same 07 range for both the 
LLA predictions and the N-body simulation results. For all the negative power models, as easily 
seen from our results in Figure 1, the agreement between the Local Lagrangian predicted P(rj) 
and the N-body simulation P s (rj) is reasonably good in the range 0.35 < ai < 0.72 and it starts 
breaking down beyond that point, the Local Lagrangian P(rj) systematically overestimating the 
peak of P s (rf) and underestimating the tail. However, the picture is quite different for the n = 
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and n = +1 power spectra. Here we have agreement between P(rj) and P s (r}) only at very low o\ 
values. Our results indicate not only that the Local Lagrangian Approximation breaks down very 
early in the linear regime for these values of n, but also that its range of validity decreases as the 
power index moves toward more positive values. 

In order to examine the behavior of the higher moments of the density field, we calculated the 
quantity 

^pi, (i5> 

where the nonlinear < 5 2 > was obtained from our simulations. We observed a decrease with 
increasing n, a result in agreement with Lokas et al. (1996) and explained as the effect of 
previrialization. Furthermore, our N-body simulations indicate that a similar behavior is also 
exhibited by the skewness-related quantity 

<6 3 s >-faf , , 

6 • ( 16 ) 



4. Conclusions 



Despite its simple form, the Local Lagrangian Approximation yields the right predictions for 
the evolved P(rj) in the weakly nonlinear regime (07 < 0.72), but only for the cases where n < — 1. 
For n > 0, the LLA fails to agree with the simulations except at very early stages in the evolution. 
Our results are in rough agreement with those of Scoccimarro &; Frieman (1996), who examined 
the contribution of next-to-leading order terms in the perturbative expansion of (cr 2 ). They found 
that these terms diverge for n > — 1 and converge for n < — 1. Our results support the conclusion 
that tree-level perturbation theory fails for n > 0. However, we find that tree-level theory, as 
expressed in the LLA, can be applied in the case n = — 1. This is not totally contradicted by the 
results of Scoccimarro & Frieman (1996), because in this case the divergence of the next-to-leading 
order term in their calculations is only logarithmic. 

To the degree that P(rj) carries all the information about its higher order moments, one 
would expect the predictive power of the LLA to extend to the calculation of such moments in its 
range of validity. Calculation of the skewness based on the LLA and including one loop corrections 
(Scoccimarro 1996) leads to reasonable agreement as expected. However, the smallness-of-c) 
condition is probably not satisfied for the power spectra n = 0, +1 we examined. Possible reasons 
for the breakdown of the perturbative approach may be the strong coupling of the long k modes 
due to non-local interactions, as well as multistreaming which may alter the picture even at the 
weakly nonlinear stage (Bharadwaj 1996). 
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^From a practical point of view, the Local Lagrangian Approximation seems to provide 
an excellent prediction for the evolution of the density PDF in the range n < — 1. Since the 
power spectrum in the quasilinear regime is close to n = — 1 (Klypin and Melott 1992), the LLA 
may be useful in comparing with observations. In particular, it may be possible to invert our 
mapping to go from the evolved (observed) PDF backwards to the initial conditions. Since the 
LLA provides better agreement with the evolved PDF than the Zeldovich approximation, it is also 
worthwhile to determine whether a modification of the Zeldovich approximation can be derived 
which corresponds to the mapping in equation (||). 
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V V V 

Fig. 1. — Comparison of the Local Lagrangian Approximation top-hat smoothed density P(rj), shown as a 
continuous line, to the N-body simulation results, shown as points with 1— a error bars, for r\ = p/p < 4 and 
— 3 < n < — 1. For each power index n we show three different a^ inear -) regimes corresponding to smoothing 
lengths A. 
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Fig. 2.- 



Same as in Fig.l but for spectral indices of n = 0, +1 respectively. 



-10- 



REFERENCES 



Bernardeau, F., 1992, ApJ, 392, 1 
Bernardeau, F., 1994, A & A, 291, 697 
Bernardeau, F., & Kofman, L., 1995, ApJ, 443, 479 



Bharadwaj, S., 1996, preprint, astro-ph /9606 1 2 1 



Juszkiewicz, R., Weinberg, D.H., Amsterdamski, P., Chodorowski, M., & Bouchet, F., 1995, ApJ, 
442, 39 

Klypin, A. A. & Melott, A.L. 1992 ApJ, 399, 397 

Kofman, L., Bertschinger, E., Gelb, J.M., Nusser, A., & Dekel, A., 1994, ApJ, 420, 44 

Lokas, E.L., Juszkiewicz, R., Bouchet, F.R., & Hivon, E. 1996, ApJ, 467, 1 

Melott, A.L., & Shandarin, SF. 1993 ApJ, 410, 469 

Melott, A.L., Splinter, R.J., Shandarin, SF., & Suto, Y. 1996 preprint 

Protogeros, Z.A.M, & Scherrer, R.J. 1996, MNRAS, in press, |astro-ph/9603i~55| (PS) 

Scoccimarro, R. 1996, ApJ, submitted 



bcoccimarro, R. k Frieman, J. 1996, ApJ, submitted, |astro-ph/9602070| 



This preprint was prepared with the AAS macros v3.0. 



